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A METHOD FOR PREDICTING THE LOW-SFEED CHORDWISE PRESSURE 
DISTRIBUTION OVER SHARP-EDGE AIRFOIL SECTIONS WITH 
PLAIN FLAPS AT THE LEADING AND TRAILING EDGES 
By Robert J. Nuber and Jones F. Cahill 


SUMMARY 


A method is presented for the calculation of the low-speed chord- 
wise pressure distribution over various sharp-edge airfoils equipped 
with leading-edge and trailing-edge flaps of arbitrary size and deflec- 
tion. The method was developed from an analysis and correlation of 
experimental pres sure -distribution data obtained from an investigation 
of a sharp-edge airfoil having leading-edge and trailing-edge flaps. 

The generalization of the data obtained from tests of a single airfoil 
to apply to other relatively thin sharp-edge airfoils rests on the funda- 
mental assumption that for sharp-edge airfoils the separation phenomena 
controlling those components of the pressure distribution which can not 
be calculated from potential -flow theory do not vary appreciably with 
variations in the detailed shape of the airfoil. 


INTRODUCTION 


Thin sharp-edge wings designed to minimize wave resistance have 
been proposed for use on high-speed aircraft. As a result, leading- 
edge high-lift devices are required for satisfactory low-speed charac- 
teristics and the necessity has arisen of obtaining adequate section 
load data for such lift ' augmenters in conjunction with trailing-edge 
flaps. An investigation has accordingly been made in the Langley two- 
dimensional low-turbulence tunnel of the aerodynamic loads over leading- 
edge and trailing-edge plain flaps ( 0.15 chord and 0.20 chord, . respec- 
tively) on a 6 -percent -thick symmetrical circular-arc airfoil section 
(reference 1) . 

In an effort to provide the designer- with additional section load 
information, a generalized method has been developed from an analysis 
of data of reference 1 which permits the determination of the chordwise 
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pressure distribution over sharp-edge airfoils equipped with leading- 
and trailing-edge flaps of arbitrary size and deflection. The develop- 
ment of the method is presented in this paper together with an illustra 
tive example. 
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COEFFICIENTS AND SYMBOLS 


airfoil-section lift coefficient ( — 

\qc 


airfoil-section ideal lift coefficient; lift coefficient at 
which stagnation point occurs at leading edge 


change in ideal lift coefficient caused by flap deflection 


airfoil-section additional lift coefficient due to angle of 
attack (c* - C2 bS ) 


surface pressure coefficient 
ible flow) 



(in incompress- 


pressure-difference coefficient across airfoil 


PU ‘ pl' 


flap-section normal -force coefficient 
flap-section hinge-moment coefficient 



airfoil lift per unit span 
free -stream total pressure 
local static pressure 

flap normal force per unit span, positive upward 

flap hinge moment per unit span, positive when trailing edge 
tends to deflect downward or leading edge upward 

airfoil chord 

flap chord 
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free -stream dynamic pressure 


PT 

2 


P 

V 

dr 


V 

V' 

Avo 


free -stream density 
free-stream velocity 

airfoil-section angle of attack, degrees 

flap deflection, positive when deflected below chord line, 
degrees 

local velocity on airfoil surface 

incremental local velocity on airfoil surface due to separation 

additional local velocity on airfoil surface due to departure 
from the ideal lift coefficient 


x distance behind leading edge, inches 

- (?) ' 

R Reynolds number 

V kinematic viscosity 

Subscripts: 

N leading-edge flap 

F trailing-edge flap 

i ideal 

U upper surface 

L lower surface 

b 5 * refers to conditions' at the ideal lift coefficient with flap 

deflected 

a refers to difference between conditions at the ideal lift 

coefficient and any arbitrary lift coefficient 
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DERIVATION OF METHOD 


Velocity distributions as calculated by potential -flow methods 
generally bear little resemblance to those obtained experimentally on 
sharp-edge airfoils because of the existence of extensive regions of 
separated flow. When analysing the velocity distribution about sharp- 
edge airfoils with flaps, however, it is convenient for the resultant 
distribution to be broken down into various component parts as is done 
in the case of airfoils in potential flow (reference 2). The most gen- 
erally used breakdown considers the resultant velocity distribution to 
be made up of the following three components: 

(1) Distribution of velocity about the basic symmetrical airfoil 

at zero angle of attack v/V 

(2) Incremental additional velocity distribution due to departure 

of the airfoil from the ideal lift coefficient Av a /V 

(3) Mean-line velocity distribution 

(a) Caused by airfoil camber Av/V 

(b) Caused by flap deflection (Av/vJ^g 


In the present paper, the only type of mean-line velocity distribution 
considered is that resulting from flap deflection, since the data used 
in the analysis are for a symmetrical airfoil section. 

In terms of the three -component velocities, the complete velocity 
distribution about an airfoil at any lift coefficient is given approxi- 
mately by: 




(1) 



(2) 


For the basic thickness form at zero lift, the velocity distribution 
can, in any case, be calculated by .the methods of references 3 a nd k. 

In the absence of flow separation, the component Av a /V is usually taken 
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to be- a linear function of the additional lift coefficient c i & , that 

is, the difference between any arbitrary lift coefficient and the ideal 
lift coefficient, and can be calculated by thick-airfoil theory. If 
extensive regions of separation do not exist, the component ' Av/V 
resulting from airfoil camber or flap deflection can also be calculated. 
The methods of thin-airfoil theory (references 5 «nrt 6) are usually 
employed for this purpose. . 

For sharp-edge airfoils on which flow separation limits the appli- 
cability of potential -flow methods, the problem of developing a general 
method of determining the velocity distribution resolves itself into a 
determination of the manner in which the various component distributions 
vary with c^ and 8. First, the manner in which the additional 

velocity distribution Av a /V varies with the additional lift coefficient 

must be found-. The use of the experimental pressure-distribution data 
of reference 1 and the following relation provides the solution: 


A v a _ \[S\j - \/Sl /Av\ 

v 2 ■ w ; b5 


Next, the extent to which > the theoretical velocity distribution due to 
flap deflection is realized experimentally must be determined. ' In order 
to determine the variation of Av a /V with lift coefficient and the com- 
parison of the experimental and theoretical velocity distributions, the 
ideal lift coefficient must be known. For any combination of leading- 
edge and trailing-edge flap deflections, the ideal lift coefficient can 
be calculated by the methods of reference 5; however, because of flow 
separation, a correlation must be made between the theoretical and 
experimental ideal lift coefficients. Finally, it becomes necessary to 
determine the velocity distribution about the basic thickness form at 
zero lift, that is, by definition: 



The value of v/v for the symmetrical airfoil at zero lift can, of 
course, be calculated by potential— flow methods^ however, the extent of 
the separated flow on the upper and lower surfaces and, therefore, the 
effective value of v/V which must be used in equations (l) and (2) 
vary with lift coefficient and flap deflection. Consequently, the value 
of v/V - for the symmetrical airfoil at zero lift, determined theoreti- 
cally, must be corrected by an increment v'/V which" is a function of 


I 
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lift coefficient and flap deflection. The value of v'/V will be 
determined from the data of reference 1. 

Although the method to be presented is developed on the basis of 
data from tests of only one airfoil section, it is thought that the 
results may be generally applied to other relatively thin sharp-edge 
airfoil sections because, for such airfoils, the separation phenomena 
controlling those components of the pressure distribution which cannot 
be calculated from potential- flow theory, would not be expected to vary 
appreciably with variations in the shape of the airfoil. 

Ideal, lift coefficient .- As indicated previously, a knowledge ,of 
the association between the theoretical and experimental ideal lift 
coefficients resulting from flap deflection, c Zgg, is essential. The 

change in ideal lift coefficient is equal to the sum of two-component 
changes, one resulting from leading-edge -flap deflection and the other 
from trailing-edge-flap deflection. Each of these components may be 
calculated separately and added linearly. For each leading-edge- and 
trailing-edge-flap deflection investigated in reference 1, the lift 
coefficient at which the stagnation point occurred at the leading edge, 
c^g, has been determined, from the experimental data. The results are 

compared in figure 1 with those calculated from thin-airfoil theory. As 
shown in figure 1 the experimental values of c'j vary almost linearly 

with trailing-edge-flap deflection above 10° and, for the leading-edge 
flap, the theoretical coefficients are identical with those obtained 
experimentally. In calculating the ideal lift coefficients cygg^ 

resulting from deflection of a leading-edge flap, the theoretical value 
may, therefore, be used. 

For trailing-edge-flap deflections above 10°, the experimentally 
determined values of the ideal lift coefficient cjggp, are considerably 

lower than indicated by the theory. In order to provide a means of 
determining the change in ideal lift coefficient associated with deflec- 
tion of trailing -edge flaps of different chords, Allen’s equation, pre- 
sented in reference 7, was used. In this equation the ideal normal- 
force coefficient is expressed as a function of flap-chord ratio, flap 
deflection, and pitching -moment increment resulting from flap deflection. 
Values of c^gg^ (the normal-force coefficient was taken to be essen- 
tially the same as the lift coefficient) obtained by this method are 
plotted in figure 2 against trailing-edge-flap deflection for flap-chord 
ratios ranging between 10 and 50 percent. The values of the quarter- 
chord pitching-moment increment required for the determination of. these 
curves were determined from numerous experimental data. -These ideal 
lift coefficients (fig. 2) represent average values obtained from a 
series of tests of plain flaps on a large number of conventional airfoil 



NACA RM L50H17a 


7 


sections. Similar computations were also made for the 20-percent-chord 
flap on the circular -arc airfoil used. in the present analysis and, as 
expected, the results were identical with the corresponding data in 
figure 2. For any profile with plain flaps, therefore, the results of 
figure 2 can he used for the determination of the ideal lift coefficient. 

Mean-line velocity distribution .- The distribution of velocity 
resulting from flap deflection (Av/V )h§ was computed from the data of 
reference 1 for various flap deflections by means of the following 
equation: 





This expression was obtained by subtracting equation (2) from equation (1) 
since, by definition, Av a /V is zero at the ideal lift coefficient. The 
data thus obtained for various deflections of the leading-edge flap were 
found to be very nearly independent of flap deflection when expressed in 

^ / c Zb6u‘ A comparison of the mean value of 


the form of 


(v) b 6 N / C 


V /b6]y 


c 2b&jj plotted against percent chord as determined by theory and 


experiment (fig. 3(a)) shows good agreement. It is concluded, therefore, 
that the mean-line velocity distribution resulting from deflection of 
leading-edge flaps of various chords can be calculated theoretically 
with a sufficiently high degree of accuracy. Due to the effects of 
separation near the trailing edge, however, the experimental velocity 
distributions resulting from deflection of the plain trailing-edge flap 
differed markedly from those predicted by the theory, particularly for 
large flap deflections. A different distribution for each trailing-edge - 
flap deflection (reference 1) was determined, therefore, and the results 
are presented in figure 3(b) in the form of (Av/V)- b5 ' against percent 


chord . 


As a basis for extending the analysis to apply to sharp-edge air- 
foils having trailing-edge flap-chord ratios other than 0.20, the 
normal-force distribution was determined from the pressure 

distribution at the ideal lift coefficient for several trailing-edge- 
flap deflections by the following relation 
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cz bS F . CZ bSp CZ b&F 


When compared with the distributions presented in reference 7 for a 
0.20c trailing-edge flap, good agreement was obtained. It is probable, 
therefore, that the normal-force distribution Pb&p/c^g^, and conse- 
quently the velocity distribution, may be determined with satisfactory 
precision for any desired trailing-edge flap-chord ratio and deflection 
from table III of reference 7 . 

Additional velocity distribution .- The values of the local incre- 
mental velocity ratio Av a /V were determined from the data of refer- 
ence 1 and equation (3) and, when plotted as a function of the addi- 
tional lift coefficient (cj a = c i - c 2^,5), were found to be essentially 

independent of either leading-edge- or trailing-edge -flap deflection. 
Average values of Av a /V are plotted against c i & in figure 4 for 

various chordwise positions. It is thought that these values of Av a /V 
(fig. 4) can be used for various flap-chord ratios since, after the 
leading edge has caused separation of the flow, any differences in air- 
foil contour behind that point would have only secondary effects 
on Av a /V. . 


Effective basic velocity distribution .- The velocity distribu- 
tion v/V about the basic symmetrical airfoil (reference 1) at zero 
lift is presented in figure 5. As previously indicated, such velocity 
distributions can be calculated by the available theoretical methods 
given in references 3 and 4. The effective values of v/V which must 
be employed in equations (l) and (2), however, vary with both additional- 
lift coefficient and trailing-edge -flap deflection because of separation 
phenomena. The increment v'/V which must be added to the theoretical 
basic velocity distribution was determined from the following relation: 



(5) 


Since v'/V is a function of both trailing-edge -flap deflection anH 


lift coefficient, it may be broken down into two components; namely, 

( v'/V )p and (v'/v) a , respectively. The values of (v'/V) F were 

determined first from equation (5) by using the experimental pressure 
distributions at the design lift coefficient. Values of the total change 

flL' 

W 


in basic velocity distribution (— ] = 


V CLJLUC 

PH 


were determined 
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from equation (5) by using the experimental pressure distributions at 
various lift coefficients. From the results thus obtained, the values 
of (v*/V)p were subtracted to obtain (v f / v ) a - It should be pointed 

out that deflection of the nose flap had no appreciable effect on the 
shape of these velocity distributions when expressed as a function 
of cz a . 


For various .chordwise positions, values of (v , /V)p are presented 
in figure 6 as a function of trailing-edge -flap deflection. Forward of 
the 40-percent-chord station values of this component of velocity were 
found to- be negligibly small. The chordwise position of (v ! /V)f is 


expressed in terms of 


_x/c 


and 


i - * 
c 


for points ahead of and 


1 - Ep “““ Ef 

behind the hinge, respectively. In this form, the results' are correlated 
so that they may be applied to sharp -edge airfoils having trailing-edge 
flaps of varying chord. This method of correlation is thought to be 
justified since the distribution of (v'/v)f is a result of separation 
at the flap hinge and has been shown (reference 7) to be similar for 
various hinge locations. 


The results of the determinations of (v'/V) a are shown in fig- 
ure 7 • As would be expected, the values are independent of flap deflec- 
tion when expressed in terms of c i & . 


APPLICATION OF METHOD 


In order to demonstrate the application of .the method derived 
above, the computations required to determine the pressure coefficient 
at 55 percent of the chord of a 6-percent-thick symmetrical circular- 
arc airfoil section with a 10-percent-chord leading-edge flap 
deflected 30° and a 30-percent -chord trailing-edge flap deflected L0° 
are worked out in the following example. The section lift coefficient 
is assumed to be I.65. 

I 

Lift-coefficient increments: 


From reference. 'J , 


c l 


b5 N 


2 sin 0 O ^N =0.628 


From figure 2, 

c *bS F = °*T72 
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Therefore 

ci = 1.65 - ( 0.628 + 0.772) = 0.25 
Airfoil basic velocity: 

From figure 5 (o r from theoretical calculations) 

^ = 1.078 

Incremental additional velocity: 

From figure 4 for ci = 0.250, 

cL 


Av a 

V 


0.034 


Mean-line velocity: 

The local incremental basic velocities are computed from reference 7 
in terms of the pres sure -difference coefficient Pgg 


From the theoretical relation (equation A-19, reference 7)> 

p b5 N = 0.4l6 

From table III (d) of reference 7 > . - 

Pb& F = '1.639 


Incremental velocity due to separation: 
For x/c = 0.55> 

r^I = 0-786 

And from figure 6, 


From figure 7 


(f), ■ 

ea 


= - 0.006 
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Pressure coefficient: 

Rewriting equation 6 of. reference 8, 



where v/V is the total of all the other v/V terms 
(excepting (Av/V ) b5 j 

and P b5 is the total of P b g N and P b&F , 

Su = 2.361 

S L = 0.202 

ACCURACY AND LIMITATIONS OF METHOD • 


In order to justify the method of correlation employed in the 
development of the present method, the calculated pressure distributions 
over the 6-percent -thick symmetrical circular-arc airfoil section and 
the integrated flap normal-force and hinge-moment coefficients for 
several individual and combined deflections of the plain leading-edge 
and trailing-edge flaps are compared with those obtained experimentally 
(reference 1) in figures 8 to 10. The flap pressure coefficients 
(fig. 8)- are plotted against the projected chordwise position of the 
flap orifices on the airfoil chord. The dispersion of the normal-force 
and hinge-moment results shown in figures 9 and 10 may be considered 
typical of the accuracy to be expected from the present method. For 
individual deflections of the leading-edge and trailing-edge flaps the 
normal force and hinge-moment characteristics, as a rule, are within 
10 percent of the experimental values. For combined deflections of the 
leading-edge and trailing-edge flaps, the predicted values of the loads 
and hinge moments over the trailing-edge flap remain within 10 percent; 
whereas, for the leading-edge flap the method tends to underestimate 
these characteristics to a larger degree, depending upon the magnitude 
of the flap deflections. 

The flap hinges were located on the lower surface of the airfoil 
and the flaps were in contact with the flap skirts so that, in effect, 
there was no leakage of air between the upper and lower surfaces. It is 
believed that changes in the vertical location of the hinge line will 
have negligible effects on the airfoil characteristics. If leakage at 



12 


NACA RM L50H17a 


the flap hinge were present, however, the effects may he such as to 
alter the separation phenomena particularly at low trailing-edge-f lap 
deflections. 

Although there may he some tendency of increased Reynolds number 
to alter the conditions of the boundary layer, the effects of scale 
. will probably be of an insignificant nature particularly in view of the 
negligible variations in section lift coefficient associated with sharp- 
edge airfoils. (See references 9 } 10, and 11.) 

CONCLUDING REMARKS 


The experimental pressure distributions obtained for a 6-percent- 
thick circular -arc airfoil with various deflections of a leading- and 
trailing-edge flap have been analysed. The components of the pressure 
distribution attributable to the basic thickness form, change in angle 
of attack, and flap deflection have been expressed in such a form that 
they may be applied to the calculation of the pressure distribution 
about various sharp-edge airfoils having different flap sizes and deflec- 
tions. The application of these results to other relatively thin sharp- * 
edge airfoils rests on the fundamental assumption that, for such airfoils 
the separation phenomena which control those portions of the pressure 
distribution that cannot be calculated by potential- flow theory do not 
vary appreciably with variations in the detailed shape of the airfoil. 
This assumption seems to be justified in large measure by the fact that 
the numerous pressure distributions obtained for the 6-percent-thick 
airfoil having various combinations of leading- and trailing-edge-flap 
deflections could' be analyzed and correlated in terms which To not 
explicitly involve the flap deflection. Comparison of pressure dis- 
tributions and integrated force coefficients obtained by the generalized 
method presented with some of the experimental data employed in the 
development of the method indicates that the methods by which the data 
were generalized give over-all results which are in reasonable agreement 
with experiment. 


Langley Aeronautical Laboratory 

.National Advisory Committee for Aeronautics 
Langley Air Force Base, Va. 
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Figure 1.- Variation of change in ideal lift coefficient with deflection 
of the 0 . 15 c plain leading-edge flap and 0.20c plain trailing-edge 
flap on the 6-percent- thick symmetrical circular-arc airfoil section 
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Figure 2.- Variation of change in ideal lift coefficient with plain 

trailing-edge-flap deflection. 
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Figure 3 -- Mean line velocity distribution for various flap deflections. 
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Figure 5* - Velocity distribution for the 6-percent- thick symmetrical 
circular-arc airfoil section, (reference 1). = 0°; 8m = = 
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Figure 8.- Continued. 
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(c) a Q = 6.1°. 


Figure 8.- Concluded. 
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(a) 0 . 15 c plain leading-edge flap. 

Figure 9.- Section flap normal- force and hinge-moment characteristics of 
a 6-percent-thick symmetrical circular-arc airfoil section for indi- 
vidual deflections of the plain leading-edge and trailing-edge flaps. 
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(b) 0.20c plain trailing-edge flap. 


Figure 9 -- Concluded. 
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0 . 15 c plain leading-edge flap. 


Figure 10.-. Section flap normal-force and hinge-moment characteristics of 
a 6- percent- thick symmetrical circular-arc airfoil section for combined 
deflections of the plain leading-edge and trailing-edge flaps. 
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(b) 0.20c plain trailing-edge flap. 

Figure 10.- Concluded. 
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